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ABSTRACT 

As part of our current programme to test ACDM predictions for dark matter (DM) 
haloes using extended kinematical observations of early- type galaxies, we present a 
dynamical analysis of the bright elliptical galaxy NGC 4374 (M84) based on ~ 450 
Planetary Nebulae (PNe) velocities from the PN. Spectrograph, along with extended 
long-slit stellar kinematics. 

This is the first such analysis of a galaxy from our survey with a radially constant 
velocity dispersion profile. We find that the spatial and kinematical distributions of 
the PNe agree with the field stars in the region of overlap. 

The velocity kurtosis is consistent with zero at almost all radii. 

We construct a series of Jeans models, fitting both velocity dispersion and kurtosis 
to help break the mass-anisotropy degeneracy. Our mass models include DM halos 
either with shallow cores or with central cusps as predicted by cosmological simulations 
- along with the novel introduction in this context of adiabatic halo contraction from 
baryon infall. 

Both classes of models confirm a very massive dark halo around NGC 4374, demon- 
strating that PN kinematics data are well able to detect such haloes when present. 
Considering the default cosmological mass model, we confirm earlier suggestions that 
bright galaxies tend to have halo concentrations higher than ACDM predictions, but 
this is found to be solved if either a Salpeter IMF or adiabatic contraction with a 
Kroupa IMF is assumed. Thus for the first time a case is found where the PN dynam- 
ics may well be consistent with a standard dark matter halo. A cored halo can also 
fit the data, and prefers a stellar mass consistent with a Salpeter IMF. The less dra- 
matic dark matter content found in lower-luminosity "ordinary" ellipticals suggests a 
bimodality in the halo properties which may be produced by divergent baryonic effects 
during their assembly histories. 

Key words: galaxies: elliptical — galaxies: kinematics and dynamics — galaxies: 
structure — galaxies: individual: NGC 4374 — dark matter — planetary nebulae: 
general 
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1 INTRODUCTION 



The standard cosmological model, the so-called ACDM 
(cold dark matter w ith a cosmological constant; see e.g. 
iHinshaw et al.l [20091 ) . has been challenged b y kinematical 
meas ure ments of dwarf an d spi ral galaxies dGentile et al.l 
20051; iMcGaugh et al] 120071 : iGilmore et ail 120071 



2007: 



Salucci et al. 



2008; but see e.g. 



Spano ct al. 2008; Kuzio dc Narav ct al 
Johansson et a l. 2009; Governato et al 



20101 ). The confrontation of the precictions of the ACDM 
with early- type galaxies (ETGs hereafter) is instead more 
uncertain. On the one hand, X- r avs (see | Paolillo et ail 



2003; 



Johnson et al 



2006; 



O'Sullivan fc Ponmanl l2004bl; Humphrey et al 

120091 ; iDas et al] 120101 ) or discrete trac- 



ers such_as_g2obular_clusters ( e.g. iRomanowskv et all 



2009; 



Woodlev et al 



Shen fc Gebhardtl 
2O10f )~ 



2010; 



20ld ; 



ISchuberth et all 
confirmed the presence of massive 
luminous systems, particularly at 
On the other hand, 



haloes in the most 
the centres of groups and clusters, 
ordinary ETGs, probed with planetary nebulae (PNe), have 
manifested discrepancies with ACDM exp ectations (see e.g. 
Romanowskv et al.ll2003l . hereafter R+03; iNapolitano et all 
20051 . hereafter N+05) which may be real or due to the 



limitations of observations and dynamical analysis. 

ETGs are diffi cult to p r obe w it h standard kinematical 
techni ques 
2004al; 



dPaolillo et al 



Pellegrini & Ciotti 



2003 



2006 



O'Sullivan fc Ponmanl 
20061 ; 



Bcrg ond et all 



Pellegrini et al.ll2007h . while they are w ithin the reach of the 



Planetary Nebula Spectrograph (PN.S: lDougfas et al.ir2002l ) 
which along with other instruments is producing large kine- 



Peng et all 120041: iMerrett et all 12006 


; 


Douglas et all 120071. 


hereafter D+07; iNoordermeer et al. 


2008; iMendez et al. 


2009: NaDolitano ct al. 2009. hereafter N+09:ICoccato et al. 


2009. hereafter C+09: iHerrmann 


& Ciardullol 2009: 


Teodorescu et al.ll2010l). 



One of the main findings emerging from these obser- 
vations is the bimodal behavior of ETC velocity dispersion 
profile s in the outer regions: s teeply falling and roughly con- 
stant (|Napolitano et alll200Sl : C+09). These profiles seem to 
generally (but not perfectly) track the bimodality of the cen- 
tral regions of ETGs, which fall into the two classes of disky, 
fast rotat ors of "ordinary" luminosity, and boxy, bright slow 
rotators (ICapaccioli et al.lll992l ; iKormendv fc Benderlll996] ; 
lEmsellem et all |2007| )~ The velocity dispersion profiles are 
shaped by the combination of orbit structure and mass dis- 
tribution, but it is still unclear which of these drives the halo 
differences between the two galaxy classes. 

In inferring the mass and the orbital structure, the dy- 
namical modelling of the PN data has been so far focused 
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sion profiles (R+03; D+07; 


De Lorenzi et al.ll2008l. hereafter 


DL+08; Dc lorcnzi et all 


20091. hereafter DL+09; N+09; 


Rodionov & Athanassoulal 


20101: cf. IWeiimans et all 120091; 


Forestell & Gebhardt 20101). N+09 summarized the results. 



comparing constraints from PNe with the ones on group 
central "bright" galaxies from X-rays and globular clusters, 
and drew the tentative conclusion that there is a strong 
transition between low- and high-concentration DM haloes. 
Such a peculiar trend could imply a transition in the role 
of baryons in shaping DM haloes, or a problem with the 
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Table 1. NGC 4374: basic data for the dynamical analysis. 



Parameter 



Value 



Reference 



R.A. (J2000) 
Dccl. (J2000) 

^sys 

(m - Af) 
A v 

My 

R c from SB fit 
R c adopted 
<J0 



12 h 25 m 03.7 s 
+12° 53' 13" 
1060 km s" 1 
31.17mag 
0.131 mag 
-22.41 ±0.10 mag 
113.5" ± 11 
72.5" ± 6 
284 km s" 1 



NED 1 

NED 

NED 

Tonrv et al. (2001) 2 

Schlegel et al. (1998) 

Sec.O 

Sec.[2~T1 

Sec.[2~T1 

HyperLeda 3 



NOTES - (1): llrttpT/neciww.i pac 



(2) 
(3) 
(4) 



corrected by -0.16 mag (sec Jensen et al. 2003) 
|http:/ /leda.univ-lyonl.fr/| jPaturel et al.l l2003) 
adopting our (m — M)q. 



ACDM paradigm itself (see also N+05, and N+09 for a de- 
tailed discussion). 

The picture is far from clear and calls for more ex- 
tensive analysis. In this paper we investigate the giant 
galaxy NGC 4374 (M84) using the stellar and PN kine- 
matics data previously presented in C+09. This is a bright 
El galaxy (~ 3L* luminosity) located in the Virgo clus- 
ter core region. It may be part of a group falling into 
the Virgo cluster, but it does not show any signs of be- 
ing a grou p-central object. Mass m odels have been con- 
structed by iKronawitter et alj (|2000| . hereafter K+00), and 
ICappellari et al. (|2006l . hereafter C+06) using stellar kine- 
matics within 1 R c (the effective radius enclosing half the 
projected light). Extensive ground-base d photometry has 
been analyzed in lKormendv et al.l (2009). 

NGC 4 374 hosts an AGN as dem onstrated by X-ray 
jet emission (iFinoguenov fc JonedlaiOll ) correlated with two 
radio lobes llLaing fc Bridlejl 19871). and connected to a mas- 
sive central black hole ( Bower et al.lll998l ) . The hot interstel- 
lar gas in the galaxy is highly disturbed and not amenable 
to a standard X-ray based mass analysis (|Finoguenov et all 

mi). 

As a representative of the "bright-ETG" population 
with a flat-dispersion profile (see e.g. C+09), NGC 4374 
provides an important opportunity to investigate the differ- 
ence between the low-concentrations inferred from PNe and 
the high-concentrations from globular clusters and X-rays. 
These tracers have so far been applied to dfferent classes of 
galaxies, which suggests the possibility that there are sys- 
tematic differences in the mass tracers themselves. Alter- 
natively, the mass inferences may turn out to be robust to 
the type of tracer used, and then should be examined in 
more detail to see if they are explainable within the ACDM 
framework. 

The paper is organized as follows. Section[2]presents the 
NGC 4374 PN system properties like radial density, velocity 
dispersion and kurtosis profile, comparing them with the 
stellar light surface brightness and kinematical profiles. We 
analyze the system's dynamics in Section [3] and discuss the 
results in relation to previous galaxy analyses in Section U 
In Section[5]we draw conclusions. An Appendix covers model 
variations with an alternative choice of rejected outlier PNe. 
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Figure 1. Distribution of line-of-sight velocities of PN candidates 
around NGC 4374, as a function of radius, and relative to the 
systemic velocity (1060 km s -1 ). Red X symbols mark objects 
designated as outliers and green boxes show the bona fide PNe. 
The dotted line shows the 3 <r velocity envelope. The dash-dotted 
line the 3 a velocity envelope corrected by the energy injected by 
the interaction with NGC 4406 (see Appendix). 



2 PN SYSTEM PROPERTIES 

The data that will be the basis of our dynamical modelling 
were presented in C+09, which can be consulted for details 
of observations and data reduction. Deep long-slit stellar 
spectra were obtained with the VLT+FORS2 spectrograph 
along the major and minor axes, and 454 PN candidate ve- 
locities with the WHT+PN.S. Observations were carried out 
on two different runs (1-4 Apr 2005, 29-3 Mar 2006) with 
quite uniform seeing conditions (~ 1.2"). To accommodate 
the anticipated kinematics range for the galaxy, filter AB 
at 0° tilt was used, which has an estimated bandpass of 
~ 5026A with 36A FWHM. 

Here we begin by revisiting some of the data charac- 
terization steps with a few differences optimized for the dy- 
namical analysis. 

We present the basic properties of the field stars and 
PNe in NGC 4374, including their distributions in space 
and velocity. Since an important assumption of our models 
is that the PNe are a fair tracer population of the field stars, 
we compare throughout the properties of the stars and PNe. 
The full line-of-sight velocity field of the PN system of the 
galaxy has been discussed in C+09 (see, e.g., their Fig. 3). 
Both the galaxy light and the PN distribution appear round 
so we will assume spherical symmetry and use, as radial 
distance from the galaxy centre, the projected intermediate 
axis R m , which is related to the semi-major axis radius R a 
and ellipticity e by Rm = R a (l — e) 1 ^ 2 [where e(R a ) is taken 
from lKormendv et al.ll2009t ]. 

For the dynamical analysis in this paper, we have con- 
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centrated on identifying possible outliers which could be 
due to unresolved background emission-line galaxies, or to 
PN pair mismatches in crowded regions or in the case of 
NGC 4374, to PNe belonging to the nearby giant elliptical 
NGC 4406. As done in D+07 and N+09, we have combined 
a 3 a cli pping criterion plus a "friendless" algorithm intro- 
duced in iMerrett et"afl (120031 ). In Fig. [1] we show the PN 
individual velocities versus R m where we have marked with 
red crosses the PNe which were either outside the 3 a ve- 
locity envelope or turned out to be friendless (i.e. having a 
velocity more than 3 a away from the average velocity of 
their 20 nearest neighbours). 

Using this approach we exclude 6 out of 457 PNe from 
the C+09 catalog. Some of these 6 differ from the outliers 
identified by C+09 because the friendless algorithm is now 
applied to the raw data set rather than a point-symmetrized 
version. The outliers show a notable asymmetry with respect 
to the systemic velocity, which motivates the use of the non- 
symmetrized friendless algorithm, and which is probably due 
to a fly-by encounter of a nearby giant galaxy as discussed 
in the Appendix. 

The outlier selection is not foolproof and is a potential 
source of bias in the analysis. 

In the Appendix, we also explore the impact on the dy- 
namical models of varying the outlier selection, and find that 
the mass results are not significantly affected by changes 
in the outlier selection, while the anisotropy inferences are 
sensitive to the classification of a small number of objects. 
Follow-up spectroscopy of these objects would clearly be 
valuable. 

We next examine the spatial distribution of the final 
catalog in £|2.1l and the velocity dispersion and the kurtosis 
in $TZ\ 



2.1 Surface photometry and PN spatial 
distribution 

For t he galaxy light, we hav e used the surface photometry 
from iKormendv et al' (2009) as in C+09, but we have re- 
duced the major/minor axis to a single profile as a function 
of R m , as shown in Fig. [5] [hereafter, R m and R will be used 
interchangeably for the intermediate-axis radius]. 

To characterize the stellar luminosity profile, we 
parametrize the surface brightness (SB) profile by the Sersic 
law: 



fi(R)-n(0) oc (R/a s ) 1/r ' 



(1) 



where as is a scale length and m describes the "curva- 
ture" of the profile (|Sersidll968l ), Fitting the dataset from 
R = 1" to 465", we find as = 0.00003", m = 6.11, and fi = 
9.74 mag arcsec -2 , and R c = 113.5". Other values obtained 
the literature include 51" fr om | de Vaucouleurs et ail 



(|l99lf ). 53" from iBlakeslee et alj 
93" from 



(2001 



72.5" from C+06, 
and 142" from 

(jKormendv et all 12009]; along th e semi-major axis), and 



Hopkins et all (|2009h , 114 



204" from lJanowiecki et alj l^Old ). 

These differences do not mean that the galaxy's lumi- 
nosity profile is not reasonably well known over the region 
which we will be modelling, but that for high-Sersic-index 
galaxies, certain characteristic quantities such as i? e and 
total luminosity require considerable extrapolation and are 
poorly constrained. 



This is not a problem that we will solve overnight, and 
for the sake of using an R c parameter that is equivalent 
to the most common usage in observations and theory, we 
adopt R c = 72.5" from the wide-field R 1 ^ 4 growth-curve fit- 
ting of C+06. This differs from our approach in C+09, where 
for th e sake of uniformity we adopted the IBlakeslee et al.l 
|200ll ) values, which will not be reliable for very extended 
galaxies because of the narrow imaging fields used. Our mod- 
elling will all be conducted in physical units, so this choice of 
R e impacts only the quoting of radial ranges in some cases, 
and the comparisons to simulated galaxies. 

With our full Sersic solution, the extinction-corrected 
total luminosity in the l/-band is 7.64 x 10 10 Lv,q, or My = 
—22.4; the uncertainties in the outer surface brightness 
profile yield a (model-dependent) total luminosity uncer- 
tainty of ~10-15%. These and other global parameters for 
NGC 4374 are listed in Table [T] For practical use in the 
Jeans modeling, we have also produced a smoothed density 
profile from the data made by a combination of a simple 
interpolation of the data up to 290" and our Sersic model 
outside this radius. 

We next compare the spatial density of the PNe with 
the field stars, using the PN number density complete to 
m* + 1.1 (see C+09, Table 7). Note that while C+09 used 
R a , we bin the data using R m . 

Given an arbitrary normalization, the PN profile 
matches the stellar photometry remarkably well (Fig. [2]) — as 
also generally found in a larger sample of galaxies by C+09. 



2.2 The dispersion and kurtosis profiles 

The rotation and velocity dispersion along the major and 
minor axes of NGC 4374 have been discussed in C+09 (their 
fig. 7), together with the 2D radial velocity field (their fig. 3). 
For the spherical analysis in this paper, we reduce these data 
to a single average velocity dispersion profile after having re- 
scaled the two axes to the intermediate-axis radius R m . 

To obtain the azimuthally averaged profile, the rotation 
and true dispersion profile are folded into an root-mean- 
square velocity profile vbms = where v and a 
are the rotation and dispersion components respective!^. 
This RMS velocity is a measure of the total kinetic energy, 
and we henceforth loosely refer to it as the velocity disper- 
sion or VD. We combine the stellar data from the differ- 
ent axes by averaging, while folding the (small) systematic 



1 In the following we will use spherical Jeans equations for non- 
rotating systems. Although NGC 4374 has no significant rotation, 
the use of the i>rms wU 1 ensure that there is no rotation contri- 
bution missing in the equilibrium balance. 

2 In the long-slit stellar data, v and a are not the true clas- 
sical moments but fit parameters in a Gauss-Hermite series 
which includes the higher-order moments /13 and /14. In prin- 
ciple, we should convert these fit parameters into revised es- 
timates of the classical mom ents, e.g. using equation 18 of 
van der Marel fc Framd dl993h . Doing so would lower the outer 
stellar dispersion profiles by ~ 10%. However, it is notoriously dif- 
ficult toextra£treliable measurements of higher-order moments 
(e.g. IShapiro et al]|2006h . and we are not confident that the /14 
measurements in this case are accurate. To avoid introducing spu- 
rious corrections to the kinematics, we therefore assume the v and 
cr fit parameters are good estimates of the classical moments. 
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Figure 2. Radial surface density profiles of the field stars (V- 
band; blue star symbols) and of the PNe (green squares) in 
NGC 4374. The PN number counts have been corrected for spa- 
tial incompleteness, and arbitrarily normalized to match the stel- 
lar data. The vertical error bars of the PN data in this and in the 
following figures represent the 1 a uncertainties (based in this case 
on counting statistics and completeness correction uncertainties), 
while the horizontal error bars show 68% of the radial range of 
the PNe in each bin. The purple curve is a Sersic model fit to the 
stellar photometry, and the gray solid curve is the interpolating 
profile. The vertical dashed lines show the spatial completeness 
interval of the PN system. 

differences into the final uncertaintiefl The PN VD is cal- 
culated using a classical expression for the variance of the 
discrete velocities around the systemic velocity. Note that 
the rotation amplitude of ~ 50 km s _1 is not dynamically 
significant compared to the dispersion of ~200-250 km s _1 . 

The resulting "dispersion" data are plotted in Fig. [3] 
Overall, the use of the full PN sample in the azimuthal 
averaged profile allows us to map the kinematics of the 
galaxy out to ~ 340" which is 20% farther out than the 
major/minor axis analysis performed in C+09. 

The dispersion decreases sharply from the centre out 
to 50" where the VD from the long slit data flattens at 
~ 220 km s" 1 . The PN data are consistent with the stellar 
absorption estimates in the region of overlap and possibly 
show a rise of the VD profile from 100" with a peak of 
~ 240 km s _1 at 170" (corresponding to about 15 kpc for our 

3 The uncertainties in the PN dispersion use a classical ana- 
lytic formula that assumes a Gaussian distribution, i.e. Avrms ~ 
yJHiV? /27V 2 . We expect th is approximation to pro duce accurate 
results in realistic systems (Napolitano et al. 2001)), and we have 
carried out additional Monte Carlo simulations of a simplified 
galaxy with radial orbits, finding that the dispersion is very ac- 
curately recovered with our estimator, with a possible bias to be 
~ 5% too high. 



adopted distance) and a subsequent decrease to the original 
value of 220 km s _1 where the VD stays fiat out to the last 
data point (~ 340" or 27 kpc, i.e. ~5 7? c ). This makes NGC 
4374 a prototypical system with a flat dispersion profile, 
although the uncertainty on the R a estimate (e.g. for R c — 
204" the last data point is at ~1.8i? c ) provides a warning 
that we may not be sampling far enough from the center to 
probe the full dynamical range of the system. More extended 
data (ideally in the direction opposite to the nearby galaxy 
NGC 4406 where it is more likely that the stellar kinematics 
is undisturbed) would clarify whether the velocity dispersion 
remains flat farther out, starts to decrease as observed in 
the intermediate luminosity sample or begins to rise as the 
cluster potential is probed. 

The vrms shows a bump at around 15" which is also 
seen in the major and minor profiles and might be related 
to some kinematical substructurtQ which is not evident in 
the photometric profile (see Fig. |2J. As we are mainly in- 
terested in modelling the galaxy outskirts, the presence of 
these wiggles in the kinematical profile will not affect our 
analysis. 

The nearly flat dispersion profile in Fig. [3] corresponds 
to an asymptotic slope of —0.07 ± 0.07 which is in clear 
contrast with the decreasing profiles found in intermedi- 
ate luminosity galaxies, with typical po wer-law exponents 
of -0 .2 to -0.6 (see R+03; D+07; N+09). lNapolitano et ail 
2008 and C+09 identified a possible dichotomy of early-type 
galaxies based on these dispersion slope differences, and we 
will here investigate further the dynamical implications for 
NGC 4374. 

We next consider higher-order velocity information. We 
quantify the shapes of the stellar and PN line-of-sight veloc- 
ity distirbutions (LOSVDs) in NGC 4374 us ing a classical 
dime nsionless kurtosis, k = i> 4 /(« 2 ) 2 — 3 (see ljoanes fc Gflll 
Il998l for exact expressions and uncertainties ■ Broadly 
speaking, we can expect that « ~ is a fair indication of 
isotropic orbits, « < is pertinent to tangential orbits and 
k > for radial orbits. 

In Fig. |3] we have combined the PN estimates with 
the stellar equivalent by co nverting the long slit stellar 
Gauss-Hermite coefficient Ha (Ivan der Marel fc Franxlll993l ; 
lGerhardll 993) into kurtosis estimates using the approximate 
relation k ~ 8%/67i4. 

The PN kurtosis is consistent with the stellar proper- 
ties in the region of overlap. Thanks to the large statistical 
sample, the PN data points show error bars which are fairly 
similar to stellar estimates, based on the best quality stellar 
absorption line data. The total kurtosis profile is consistent 
with zero at all radii and has a median (calculated over all 
datapoints) of 0.05 ± 0.19. 

Our previous analyses of NGC 3379 and NGC 4494 indi- 
cated global k ~ +0.2 and +0.6, respectively. However, most 
of this difference is driven by the data inside R c , where pre- 
vious work with larger galaxy samples has indicated that any 
correlations between the fourth moment and other galaxy 

4 This may be related to th e central dust rin g clearly seen in 
optical imaging of the galaxy Jjaffe et al.lll994). 

5 Monte Carlo simulations based on iNapolitano et al.l teOOlT ) 
models have demonstrated accurate recovery of the kurtosis using 
our estimator, with a systematic deviation of no more than ~ 0.1, 
see also N+09. 
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Figure 3. Composite projected RMS velocity and kurtosis pro- 
files of NGC 4374, with data from stars (filled star symbols) and 
PNe (open circles). Separated profiles of rotation and true dis- 
persion can be seen in C+09. 



prop erties are subtle (|Bender et al.l 1 19941 ; iKrainovic et al.l 
l2008h . 

In the outer parts, all three galaxies are similarly con- 
sistent with zero kurtosis, and it will be interesting to see if 
any patterns emerge with a large sample. However, as we will 
see in the next Section, interpreting the orbital anisotropy 
implications of the kurtosis requires detailed modelling. 



3 DYNAMICAL MODELS 

We present a suite of Jeans dynamical models following the 
same scheme as in N+09, to which we refer the reader for 
more details of the analysis. We will combine the photomet- 
ric and kinematical data for the stars and PNe in NGC 4374 
into integrated models in order to derive the mass profile 
and the orbital distribution of the galaxy and finally test 
whether or not it hosts a massive dark halo compatible with 
the ACDM predictions. 

Although there are other dynamical procedures such 
as Schwarzschild's m ethod and made-to-m easure particle 
methods (e.g. R+03; IChaname et all 120081 ; DL+08) that 
have been applied to discrete velocity data and are more ro- 
bust than our Jeans approach, the latter is computationally 
faster and somewhat more intuitive. Furthermore it allows 
a much larger flexibility on the range of galaxy potentials to 
be used. In the following we briefly remind the main steps 
of our dynamical procedures. 
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Figure 4. Composite projected velocity dispersion profile of 
NGC 4374, with data from stars (filled star symbols) and PNe 
(open circles). The black solid curve shows the pseudo-inversion 
mass model to fit the PN data outside 10" for the isotropic case, 
with the shaded regions showing the 1 a significance of the fit. The 
short dashed blue curve shows the solution for /3 = 0.5, the dot- 
dashed green curve the one for ft = —0.5. The long-short-dashed 
violet line shows the solution for the cosmological motivated 0(r) 
profile as in Eq. [5] The thick red solid line shows the heuristic 
/3(r) model adopted in Sect. [5~5l 



In the different formulations of the Jeans equations we 
will assume spherical symmetry. This is a reasonable ap- 
proximation because the round and boxy stellar isophotes of 
NGC 4374 (average ellipticity (e) = 0.13 and (q 4 ) = -0.4; 
see C +09), and the small V/a — 0.03 (|Cappellari et all 
2007jf| make the system a typical boxy-slow rotator which 
is highly unlikely to be very flattened intrinsically. 

Another basic assumption of our analysis is that stars 
and PNe are all the drawn from the same underlying dy- 
namical tracer population, which is well motivated by the 
agreement between the stellar and PN properties ( £|2. H and 
1 2 . 2 p . We will also in general omit the stellar kinematics data 
inside 10" from our model fits, since there appears to be 
a strong dynamical change in the nuclear region which our 
smooth Jeans models are not designe d to reproduce (par - 
tially produced by a massive black hole; iBower et ail lTggd 7 !) . 



We begin with a simple non-parametric model in g3TJ 
then introduce multi-component mass-models in H3.2I and 
additional dynamical methods in H3.3I The multi-component 
results are presented in i]3.4H3.6l and the mass profiles sum- 
marized in H3.7I 



6 Their fig. 3 illustrates an estimated family of deprojections for 
this ga laxy, with the most flattened solution having e ~ 0.2. 

7 Here lBower et ah! l ll99Sl) estimate a black hole mass of Mgjj ~ 
1.5 X 10 9 Afo which implies a sphere of influence of radius ~ 
1.7", where we have defined as the radius where M*(r < r%) = 
2Mbh, with M*(r) corresponding to the Kroupa IMF. 
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3.1 Pseudo-inversion mass model 

We start with a phenomenological approach introduced in 
R+03 and followed in D+07 and N+09, used to convert the 
observed kinematics into a mass profile M(r). This approach 
has the advantage that it is computationally light, does not 
involve Abel inversion integrals, and does not assume any 
form for M(r), nor a stellar M/L value (which will be dis- 
cussed later in this Section). A disadvantage is that it does 
not allow a direct test of any theoretical prediction (which 
we will do in the next Sections). 

For the benefit of readers not familiar with this proce- 
dure, we summarize in the following its basic steps: 

(i) Adopt a simple smooth parametric function for the 
intrinsic radial velocity dispersion profile: 



,(r) 



CO 



1 + 



(2) 



> r + ro , 

where <to, ro, r\ are a minimalistic set of free parameters. This 
model is adopted to reproduce the flat dispersion profile in 
the outer galaxy regions and is different from those adopted 
in D+07 and N+09 which were constructed to match steeply 
decreasing velocity dispersion profiles. 

(ii) Assume a given anisotropy profile, often constant or 
parametrized as a simple function: 



P{r) = 1 - a e /a r , 



(3) 



where <jg and ay are the spherically-symmetric tangential 
and radial components of the velocity dispersion ellipsoid, 
expressed in spherical coordinate^. 

(iii) Project the line-of-sight components of the 3-D ve- 
locity dispersions oy and erg for comparison with the line-of- 
sight velocity dispersion data o\ os (R). 

(iv) Iteratively adjust the free parameters in Eq. [5] to 
best fit the model to the observed dispersion profile. 

(v) Use the best-fit model (Eq. [2]) in the Jeans equation 
4-55 of lBinnev fc Tremaine! l| 19871 ') to calculate M(r): 



MM = 



G 



dlnr 



dlnr 



(4) 



where j*(r) is the spatial density of the PNe, and corre- 
sponds to an Abel deprojection of a smoothed density law 
as in H2.ll Additional quantities may then be computed, 
such as the cumulative M/L. 

Starting with the isotropic case (/3 = 0), we find that 
the simple model ([2]) is able to fit the dispersion data well 
(Fig. with some systematic discrepancies at ~ 40" that 
we will improve upon with more complicated models be- 
low. The resulting M/L profile increases steeply with the 
radius (Fig. [S]), providing a strong indication for the pres- 
ence of an extended DM halo. Note that the shaded regions 
in Figs. [4] and [5] along with the various uncertainties quoted 
below account for the 1-cr statistical confidence region in the 
parameter space (ao,ro,rj) of the dynamical model. 

The central dynamical (M/L)v = 6.5 can also be 
compared with independent stellar population anal yses of 
the stellar M / L, T «. Assuming a iKroupal (|200ll ) IMF, 
iTortora et ail (|2009h found T. ~ 3-4.5 T ,y, while 



8 Due to the modest rotation of the galaxy, we expect the spher- 
ical approximation not to cause any significant systematic issues. 
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Figure 5. Cumulative V*-band mass-to-light ratio (M/L) of 
NGC 4374 (note that the vertical axis starts from Tq y = 3). 
The curves based on the pseudo-mass inversion method are colour 
coded as in Fig. [4] We also add some of the models from the Jeans 
analysis in ^3.31 the dotted red curve is the "NFW + /3(r)", the 
dashed red curve is the same model with adiabatic contraction 
["NFW + AC + /3(r)"], and the dashed gray line is the logarith- 
mic potential model with 0(r) (see £|3,5l and ^13. 61) . The horizontal 
blue shaded region shows the stellar M/L and its uncertainty for 
the Kroupa IMF, while the green one is for the Salpeter IMF. 
The small purple shaded region is the dynamical M/L estimate 
from C+06. See text for details. 



iGerhard et all (|200l[ ) found T* ~ 4.5-6.0 T©,v (where 
we have in both cases converted from B- to V-band). 
C+06 found T„ = 3.08Tq in /-band which we convert to 
T» ~ 5.14 Tq,v after detailed comparison of the SB profiles. 
(Note that their Schwarzschild modelling analysis implies a 
dynamical T ~ 7.3 ± 0.4X0^ in the central regions, which 
agrees with our Jeans results, as shown in Fig. 0) 

We can reasonably assume T, ~ 4-6 Tq,v for a Kroupa 
IMF, which corresponds to ~ 6.5-9.5 T 0i y- for a ISalpeteil 
(|l955h IMF (see Fig. 0}. Therefore the dynamical M/L is 
suggestive of some dark matter inside R c (72.5") for the 
case of Kroupa but not Salpeter. In the following we will 
consider the stellar M/L based on the Kroupa IMF as the 
reference results, since there are arguments to consider this 
one as a universal IMF (Kroupa 2001). 

Our last datapoint (~ 340") is close to ~ 5 R c , which is 
a benchmark distance for the mass profiles (see R+03, D+07 
and N+09): here we find that the V-band M/L within this 
radius is Ys,y ~ 20 + 2 Tq,v0. The anisotropy is accounted 
for in step (ii) of the procedure by adopting constant values 

9 Hereafter we are deliberately neglecting the uncertainty on _R e 
which we have seen are unreasonably large and scale all the results 
for our assumed R c . 
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of /3 — ± 0.5 as a plausible (though not exhaustive) range of 
the stellar anisotropy. The fits to the data are just as good 
as for the isotropic case as shown in Fig. [4] 

In Fig. [5] we show the M/L(r) profiles corresponding to 
the three P values. Assuming /3 = +0.5 implies a smaller 
central M/L (~ 5.5 Tq,v) but a steeper M/L profile than 
the case of /3 = 0, while /3 — —0.5 implies a larger central 
M/L (~ 7T0,v) and a shallower M/L profile outside 1 R e 
(with the M/L consistent with the isotropic profile at all 
radii in either cases). In all cases a constant M/L is excluded 
at more than 3 a and DM starts dominating already at 1 R e , 
assuming a Kroupa IMF, and at ~ 2 R c for the Salpeter IMF 
case. 

Our outer M/L results are relatively insensitive to the 
anisotropy assumed because of a geometrical effect in certain 
regimes in radius that causes anisotropy differences to cancel 
out when proj ected t o line-of-sight velocit y dispersions (cf. 
Gerhard] 11991 F ig. 8;| van der Marellll994 Figs. 10 and 11; 



Wolf et all [2010 ). This "pinch point" occurs where the 3D 



log slopes of the tracer density profile a and the velocity 
dispersion 7 add up to (a + 7) ~ —3 (see Dek el et al. I l2005l 
Eq. 2). In a bright elliptical galaxy like N4374, the high 
Sersic index n, the large scale-length, and the flat dispersion 
profile combine to push the pinch point to fairly large radii: 
~ 100" in this case. This robustness of the mass inference 
contrasts with the case of galaxies with steeply declining 
dispersion profiles, where the mass-anisotropy degeneracy is 
particularly severe (DL+08; DL+09). 

We have also tested the anisotropy profile based on the- 
oretical expectations from merging collisionless systems as 
derived from ML05: 



P(r) = fa 



r + r a 



(•5) 



where Po ~ 0.5 and r a ~ 1.4 R c (based on the merger sim- 
ulations of D+05). Adopting this profile with r a = 101", 
we find that the VD profile matches slightly better the cen- 
tral regions, but it fits poorly the large radii datapoints. 
In this respect Eq. [S] seems to be ineffective in reproducing 
the intrinsic anisotropy of the galaxy (given the limits of 
the simple parametrization assumed in equation [2)F°I How- 
ever, the fact that radial anisotropy produces a better fit to 
the central VD while P — matches the outer parts of the 
galaxy suggests that a more complicated p(r) profile than 
the one in Eq.[S] should be applied to NGC 4374. 
For instance, looking at the kurtosis profile in Fig. [3l one 
suspects that a /3(r) profile which is isotropic in the very 
central regions (R < 5") and in the outer parts (R > 70") 
and radially anisotropic in between (R ~ 5" — 70") might 
do a better job. Following this heuristic approach we adopt 
the following formula: 



P(0 = Po 



1/2 



(6) 



where Po = 0.6 and £ = r/r a with r a ~ 30" (see also Sec- 
tion [33]). This P(r) profile is significantly different from the 
simulation-based Eq. [5] but similar to the /3(r) found from 
the detailed dynamical models of K+00 for NGC 4374 (see 

10 We tried out a wider range of r a : for smaller r a the predicted 
dispersion was still lower than the data, and for larger r a the 
dispersion progressively approached the isotropic case. 
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Figure 6. The heuristic 0{r) profile in Eq. [6] (solid line) is com- 
pared with the simulation based /3(r) ( dashed line) from ML05 
as in Eq. [5] and the modelled P(r) from Kronawittc r et aTll200ol 
(shaded region). Here only the radial range covered by the Kro- 
nawitter et al. model is shown: the matching with the heuristic 
/3(r) is good, while the ML05 formula predicts radial anisotropy 
at much larger distances from the centre. The anisotropy value 
derived from direct kurtosis inferences (see i|3.3l l is also shown 
with 1 a error bars. 



Fig. |SJ) as well as for some other galaxies in their sample 
(e.g. NGC 4278, NGC 4472, NGC 4486, NGC 5846). In 
this case the best fit to the VD is improved as shown in 
Fig. U (red curve; we will come back to this issue in H3.5p . 
The corresponding M/L profile has a central value which 
is closer to the isotropic solution (~ 6.5Tq : v) and becomes 
slightly larger outward, finally converging to the isotropic 
case asymptotically. 

The overall plausible range for the benchmark-radius 
M/L of NGC 4374 is T 5 ,v = 18-24 T©, v (including both 
statistical uncertainties as well as the systematic anisotropy 
uncertainties, given the range of P(r) profiles that we al- 
low). This result is significantly larger than the typical 
M/L found for the intermediate-luminosity galaxy sam- 
ple studied so far with the PN.S (see e.g. R+03, D+07, 
DL+09 and N+09), but more similar to the M/L esti- 
mates foun d in bright systems using globular clusters and 
X-rays (e.g. Humphrey et al.l2 006 ; Ro manowskv et "al]|2009l : 



ISchuberth et al.ll2010l : iDas et al.ll2oiol ) 

The steep increase of the M/L with radius can be quan- 
tified through the dimensionless M/L gradient (introduced 
by N+05): 



R C AT 



(7) 



Tin AT? 

where Tin is the central dynamical M/L. For NGC 4374 we 
find V^T = 0.5-0.7, which places this galaxy among the sys- 
tems with larger V<T which are discussed in N+05 as very 
dark-matter dominated. As a comparison, for NGC 3379 and 
NGC 4494 we found V^T in the range -0.05 to 0.25. 



3.2 Multi-component models: mass profiles 

The second strategy for our dynamical analysis again uses 
a Jeans analysis but begins with parameterized mass pro- 
files and projects the predicted kinematics for comparison 
to the data. Following N+09, the inclusion of higher veloc- 
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ity moments (kurtosis) in the Jeans analysis is expected to 
alleviate the mass-anisotropy degeneracy. 

In our equations, we will adopt two-component mass 
models consisting of a luminous field star distribution plus 
a DM halo. The total gravitational potential may thus be 
expressed as $ = $„ + <E>d- The stellar gravitational po- 
tential $,(r) is derived from the stellar luminosity j* (rV\ 
combined with some assumed constant T*. 

Our mass models as described below use for the DM 
either an NFW profile Ci'i.2.1 1) or a pseudo-isothermal form 



3.2.1 NFW model 

Our reference mass models aims at testing the predictions 
from simulations of collisionless DM halo formation in a 
ACDM cosmology. In this case the DM density takes the 
approximate form of an NFW profile: 



Pd(r) 



Ps 



(r/r s )(l + r/r s ) 2 ' 



(8) 



where p s and r s are the characteristic density and scale ra- 
dius of the halo. The cumulative dark halo mass is 



Md(r) = 4np a r 3 a A(r/r 3 ), 
where 

A(x) = ln(l + as) - 



1 + x 



The potential is: 



*d(r) 



4nGp s r 3 a 



In 



(—) 

\r + r a J 



(9) 



(10) 



(11) 



where G is the gravitational constant. 

The three free parameters describing the NFW mass 
model are thus Y„, p s and r a . The halo can alternatively be 
parametrized by the virial mass and concentration, M v i r = 



47rA vir p cr itr vir /3 and c v ir 



t/r a , where the critical den- 



sity is p cr it = 1.37 x 10~ Mq pc and the virial overdensity 
value is A v i r — 100. 

The expected values for these model parameters are not 
arbitrary in ACDM. For instance, in a collisionless ACDM 
universe with WMAP5 parameters, the following mean re- 
lation is expected between mass and concentratiorF 2 !: 



r(M vir 



12 



M vir 
10 n M, 



(12) 



which has a 1 er scatter of 0.11 dex, and is v alid for z = 0, 
Q.m = 0.3, f2 A = 0.7, h = 0.7, and a s = 0.8 l|Maccio et all 
2008). For comparing with models parameterized by the 
scale radius r s and density p 3 (e.g. Eq. [SJ, we find that 
Eq. [12] is equivalent to the following relation: 



11 This is obtained by Abell inversion of the observed SB in the 
central regions and the extrapolation to infinity according to the 
Sersic model of ffTLl 

12 For sake of completeness we also report here the WMAP1 



equations (see N+09 for details): 



c vir (M 
and 

P" — (lSpc 



)~18 



( M vil \ 



-0.125 



-2/3 



M 



Ps 



0.29 



lOpc 



m q pe- 



as) 



where the scatter in p a at fixed r a is a factor of 1.3. Note 
that in N+09 we used ACDM halo predictions based on 
WMAP1 parameters, which implied ~ 30% higher concen- 
trations than WMAP5. 



3.2.2 LOG model 

Our alt ernative mass mode l cons ists of a logarithmic po- 
tential (jBinnev fc Tremaind fl987l §2.2.2) which was moti- 
vat ed by observations of spiral galaxy rotation curves (see 
e.g. IPersic et al.l 1996). The potential is: 

<Mr) = ^ln(r 2 + r 2 ), (14) 

where vo and ro are the asymptotic circular velocity and 
core radius of the halo. The corresponding DM density and 
cumulative mass profiles are respectively: 



Pd{r) = 



and 



^o 2 (3rg+r 2 ) 
47rG(r 2 + r 2 ) 2 ' 



Md(r) = 



vlr 3 



Gr 2 n +r 2 



(15) 



(16) 



The three free parameters of this "LOG" model are thus Y* , 
vo, and ro. We define a virial mass relative to the critical 
density according to the same definition as in H3.2.1I (there 
is no halo "concentration" in this context). 

Unlike the NFW halo with its cuspy r _1 density centre, 
the LOG halo has a constant-density core. At larger radii, 
the density decreases as r -2 , similar to the NFW model 
near r = r„. This model allows us to maximize the stellar 
contribution to the central mass, and to test a "minimal 
DM halo" scenario. Similar models have been successfully 
used to explain the dynamic s of other galaxies of a ll types 



(e.g. iFall fc Efstathiou 1980; Bcgcman et"afl Il99ll; K+00 



Thomas et aT1l2007i ; IWeiimans et alj200St DL+08: |Pu et al.1 
20101 ). 



pc" 



3.3 Multi-component models: dynamical methods 

Our Jeans modelling approach has been extensively devel- 
oped in N+09, to which we refer the reader for the full de- 
scription of the equations adopted. Basically, in addition to 
the usual second-order Jeans equations for the velocity dis- 
persion profile, we solve the fourth-order Jeans equations 
to constrain the LOSVD with kurtosis data and reduce the 
systematic unc ertainties linked to the unknown orbital dis- 
tribution (e.g. iMagorrian fc Ballantvnel l200ll : lLokasI |2002| ; 
lLokas fc Mamonl |2003| ~ Although the higher-order Jeans 
equations are not closed in general, one can adopt a simple 
choice for the distribution function which makes the prob- 
lem tractablj^l. This simplification is arbitrary (e.g. /3 is 



13 We restrict ourselves here to functions which can be con- 
structed from the energy-dependent distribution function by 
multiplying it by a function of angular momentum f( E,L) = 
fp(E ) L~ 2 P with P = c o nst . This is a widely - used a nsatz llHenon 
Il973l; iDeionghd Il986l; IWilkinson fc Evanj Il999t lAn fc Evan" 
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assumed to be constant with radius) and does restrict the 
generality of our results, but the model is still more general 
than an assumption of isotropy. In N+09 we demonstrated 
the utility of this approach for assessing the presence of ra- 
dial orbits in NGC 4494. 

For the sake of clarity, we report in the following the 
basic steps of our analysis (for more details, see also N+09): 

(i) Set up a multi-dimensional grid of model parameter 
space to explore, including j3 and the mass profile parame- 
ters (T t ,p s ,r s ) or (T t ,v ,r ). 

(ii) For each model grid-point, solve the second- and 
fourth-order Jeans equations. 



3.4 Multi-component model results: no-DM case 



(iii) Project the internal velocity moments to eri os and 



^los 



(iv) Compute the \ 2 statistic, defined as 



-^data 

< 2 =£ 



5pf 



(17) 



where p° hB are the observed data points (<ti os and ki os ), 
pf od the model values, and 8p° hs the uncertainties on the 
observed values, all at the radial position Ri. We fit the 
PN data outside 60" (where the spatial incompleteness 
due to the galaxy ba ckground is more severe, see also 
iNapolitano et all l200ll ) and the stellar data outside 10" 
(see ©. 

(v) Find the best fit parameters minimizing the x 2 ■ I n 
practice, we find that the VD is affected by both the mass 
and anisotropy profiles, while the kurtosis is almost entirely 
driven by the anisotropy. 



One interesting side-note is that given the assumptions 
of our Jeans formalism, we showed in N+09 (Eqs. B10-B12) 
that if a system has a constant dispersion profile, we 
can estimate its internal anisotropy /3 directly from the 
data without any need for dynamical modelling. This is 
because the line-of-sight kurtosis jc is then a simple matter 
of projection effects for a given /3 and luminosity profile. 
Therefore at a radius of ~ 170" , we estimate that NGC 4374 
has an anisotropy of (5 ~ — 0.11q'|, i.e. it is near-isotropic. 

The list of mass models we will explore in the following 
Sections includes: 1) a no-DM case or self-consistent model 
where the potential is given by the stellar mass only; 2) a 
NFW dark halo to be tested against the ACDM predictions; 
3) a core logarithmic potential. The novelty of this analysis 
with respect to N+09 and all other dynamical studies on 
individual ETGs is the inclusion of the effect of the adiabatic 
contraction of the dark halo, for both the DM halo models 
as above. 



In £|3 . 1 1 we have seen that for NGC 4374, a model with a 
constant M/L with radius is ruled out by the PN velocity 
dispersion data. However, the pure-stellar potential (p s — 
or vo = 0) is the minimal model that can be tried to fit 
the dispersion and kurtosis data, allowing us to find the 
maximum stellar content of the galaxy compatible with the 
inner data points. 

The best-fit parameters of the model with an isotropic 
velocity ellipsoid (/3 = 0) are listed in Table [2] together with 
the x 2 °f the fit. 

Given the freedom to adjust T„, the model is able to 
fit the VD in the central regions ( <, 2 _R C ) with a best- 
fit T„ = 7.5 (V band). This value is consistent with the 
SSP estimates based on the Salpeter IMF, and inconsistent 
with the Kroupa IMF predictions at more than la. We will 
come back to this issue in the next Section, and note here 
that, despite the higher Y„, the no-DM model fails to re- 
produce the data since the VD falls off too quickly in the 
outer regions (Fig. [7J blue dotted line). The gap between 
the model and the data cannot be removed even by assum- 
ing extremely negative ft (see e.g. the cyan dot-dashed line 
for ft = —3 x 10 3 ) or by adopting a shallower SB profile as 
allowed by the fit errors in H2.ll 

These Jeans models are not general enough to explore 
every dynamical solution that is physically possible, but we 
judge that the data/model differences are large enough to 
render a constant M/L model highly implausible. We will 
next proceed with models allowing for the presence of a DM 
halo to find out what halo parameters are most consistent 
with the data for the two assumed DM profiles. 



3.5 Multi-component model results: NFW model 



2006), which has the advantage of being easy to integrate even 
though it does not generalize to the case of = /3(r) for the 
fourth-order moment. 



We next consider the NFW mass model (Section l3.2.1[) based 
on ACDM expectations. We initially discuss the case with 
orbital isotropy in H3.5.1I and show that this matches the 
data fairly well except near R e (namely, 20" — 100") where 
the dispersion (kurtosis) is overestimated (underestimated) 
by the Jeans models. In H3.5.1I we explore a range of con- 
stant and radially-varying f3 profiles and conclude that a 
significant radial anisotropy is ruled out at large galacto- 
centric distances, while the j3(r) profile as in Eq. |S] provides 
the best match to the data at all radii. Finally we include in 
our model the effect of adiabatic contraction in £13.5,31 and 
find that the higher central DM fraction thereby generated 
allows the data to accommodate a smaller stellar M/L, fully 
compatible with a Kroupa IMF. 



3.5.1 The isotropic model and the stellar M/L issue 

We start by assuming isotropy, and find a best fit as shown 
in Fig. [7] (green dashed), with parameters again reported in 
Table [2] ( "NFW iso"). This solution is a fairly good match 
to the data, for both the VD and kurtosis profile, which is 
a further support for the absence of strong anisotropy in 
the stellar orbital distribution. The best-fit T„ ~ 6.5T©,v 
is lower than the no-DM case because the central regions 
contain significant amounts of DM (see although it is 
still the stellar mass that determines the main kinematical 
features inside ~ 100" ~ 1.2Re. This stellar M/L value is 
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Table 2. Summary of best-fit multi-component model parameters. 



Model 


fe 1 


T„ 2 


log M, 3 




log M vir 5 / vlr 6 


/dm, 5 


T(B C ) 8 


T B5 9 


T"(-R vir ) 10 


V^T 11 


x 2 / 






( T 0,v) 


(M Q ) 




(M ) 




( T G,v) 


( T o,v) 


( T 0,v) 




d.o.f. 12 


No-DM model 


star iso 





7.5 


11.76 




11.76 





7.5 


7.5 


7.5 


(1 


123/36 


NPW model 



13.4+°t 54+^ 0.7+°;! 8+ 2 22+ 44 350 0.47 28/45 

— 0.5 —36 —0.4 —1 —8 ' 

13.3+°'^ 51 ±34 °- 8 to4 7 li 23+ -io 286 °- 65 78/40 

13.3+°'| 53 ^32 °- 8 ioI 8 il 25 il0 294 °- 72 23 / 44 

13.1+J-| 32 ^25 °- 7 io4 8 ^2 22 ill 183 °' 59 12 / 33 

13.3+°' 4 45+ 4 ° 0.7+° + 7+] 17+" 261 0.39 31/44 

— U.Q — 25 — (J. o — 1 — 1U 

13.2+°^ 3 9^24 °- 8 io 5 10 -2 32+ -14 217 lx) 40/44 

13 - 4 to'4 66 -37 °- 7 io3 7 il 18 i6 368 °- 44 15/33 



LOG model 


Model 


ft 1 


( T 0,v) 


log A/» 3 
(M e ) 


v 13 
(kms" 1 ) 


log A/ Vir 5 r+ 4 
(Mq) (arcsec) 


/dm,5 7 


T( iJ ) 8 

( T 0,v) 


T B 5 9 

( t O,v) 


r(H vlr ) 10 

( T 0,v) 




x 2 / 

d.o.f. 12 



25+J 600 °' 57 25 / 45 

24+| 485 0.73 27/44 

24+^ 440 0.65 19/33 

19+ 41 540 0.39 21/44 

24+ 7 465 0.72 28/44 

18 +f ; 414 0.47 21/33 

— 4 ' 



NFW iso 

NFW iso2 

NFW+/3 

NFW+)?(r) 

NFW+AC+iso 

NFW+AC+/3 

NFW+AC+/3(r) 



6.4 

5.5 

0.2±0.1 5.5 

0.01±0.1 5.7 

5.7 

0.30±0.15 5.5 

0.01±0.1 5.5 



11.69 9+° 

— 5 

11.62 12+i 1 

— 6 

11.62 13 lg° 
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more consistent with a Salpeter IMF than with Kroupa (to 
be addressed further in H3.5.1I) . 

The central NFW halo parameters of p s = 
0.0030t , ; l l Jg M Q pc -3 and r a = 915" ± 200" = 76 ± 17 kpc 
(see Fig. [7] which shows the joint region of permitted values 
for r s and p 3 , marginalized over the other free parameters, 
T») correspond to a virial radius, mass, and concentration 
of r-vir = 770 ± 70 kpc, A/ Vi r = (2.5±f ?) x 10 13 M© and 
Cvir ~ 9lg . These halo parameters are comfortably compat- 
ible with WMAP5 expectations (Eqs. [12] and as well 
as WMAP1 (modulo an IMF issue that we discuss below). 
Looking carefully at the details of the DM halo solution, the 
VD (kurtosis) data within 20"-100" (1.3-2 dex) are slightly 
overestimated (underestimated) by the model, which might 
be either an indication of 1) some degree of anisotropy or of 
2) a mass excess caused by a larger DM concentration not 
accounted for in the NFW halo model. 

Before we explore these two options, we will investigate 
further the IMF issue mentioned above. 

In the NFW dark halo model solutions discussed so far, 
the best-fit T, (~ 6.4) is more comfortably consistent with 
the stellar M/L predicted by the population analysis assum- 
ing a Salpeter IMF than a Kroupa IMF (see fl3Tj) . Although 
this is not a strong argument for preferring either IMF, we 
have tried to quantify the effect of Y„ on our result. 

The high T« is mainly driven by the fit to the central 
data-points and the tendency of the \ 2 procedure to favour 
more minimal halo solutions. Since our simple Jeans models 
are not designed to reproduce detailed kinematical struc- 



ture as might be present in the central regions, we lower 
the weight of the very central VD and kurtosis data-points 
(i.e. data up to 30", ~1.5 dex) in the y/ 2 minimization. In 
this case, more centrally concentrated halo solutions can be 
made compatible with the dat J*^l. Indeed, in Fig. (thin 
purple dashed line), we report the best fit obtained for the 
isotropic assumption, where a lower stellar M/L is needed, 
T„ = 5.5, which implies a dark matter halo with p s — 
0.00491^ ^3 M Q pc~ 3 and r a = 720" ± 200" = 60 ± 17 kpc 
corresponding to a virial radius, mass and concentration of 
rvir = 720 ±30 kpc, M vi r ~ 2.1 x 10 13 M Q and c vir ~ 12 (see 
also "NFW iso2" solution in Table 0. In this case, though, 
the halo concentration is higher than predicted for WMAP5 
parameters. 

In Fig. [7] it is evident that this solution has a shal- 
low velocity dispersion profile at R < 25" ~ 1.4 dex which 
is a poor match to the data and causes the high \ 2 value 
for the fit. However, the gap can be filled either with the 
presence of some (anticipated) degree of anisotropy in the 
central regions or by a DM enhancement by an adiabatically 
contracted halo. In the following, we will explore these two 
possibilities in turn. 



14 E.g., in Fig. [5] of SI3.1I a lower central M/L is found (though 
for the = +0.5 case). 
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Figure 7. Multi-component Jeans model fits to the NGC 4374 kinematics data. The stellar data are shown by star symbols, and the PN 
data are open circles. The left panels show the projected RMS velocity profiles (top) and the projected kurtosis (bottom), the right panel 
the corresponding 1 and 2 a confidence level of the p s — r s parameters marginalized with respect to T* and r a (for the "NFW+/3(r)" 
model). The curves correspond to models as in the panel legend (except "star+tan" which is not a best-fit model). The shaded regions 
on the right show the WMAP1 (gray) and WMAP5 (blue) expected region for halo parameters. The "NFW +/3(r)" model from Fig. [7] 
is plotted here for comparison with the isotropic case and repeated in Fig. [8] See text for details. 



3.5.2 Models with orbital anisotropy 

A way to produce a modelled steeper a\ os profile, for a given 
slope of the intrinsic light density profile, j t , and velocity 
dispersion (see e.g. Eq. 3]), is with some degree of radial 
anisotropy (see, e.g., Dekel et al. 2005). 

We have started with a constant anisotropy from the 
very central regions and the best-fit solution is found to 
accommodate a gentle radial anisotropy (/?o ~ 0.2) with 
a lower stellar M/L (= 5.5 Tq : v) that now agrees with 
a Kroupa IMF. The VD and the kurtosis are at last re- 
produced well at all fitted radii (Fig. [8j red dot-dashed 
line), which is reflected in an improved x 2 value in Table [2] 

("NFW+/3o")- 

The halo concentration for this solution is fairly high, 
and just consistent with the WMAP5 expectations at the 
~ la level. 

We remark here that the constant anisotropy solution 
provides a compromise model dispersion curve among re- 
gions which might have different orbital structures. For this 
reason we decided to test also the case of a radially varying 
P(r) even though our dynamical procedure is not explicitly 
designed for this. As done in N+09, we will use the kurtosis 
data to constrain /3 in the outer regions where the anisotropy 
may be approximately constant. 



Following the approach of H3.ll we use the /3(r) as in 
Eq. [6] The best-fit model is shown in Fig. [8j (black line) 
and the parameters are reported in Table[5]["NFW-r-/3(r)"]. 
The anisotropy radius r a turned out to be very close to the 
one estimated with the pseudo-inversion procedure (r a = 
33"). The match in the central regions is remarkably good 
also for the low Y„, while in the outer regions the model 
tracks the isotropic case (see left panel of Fig. [7| for a direct 
comparison), and the halo concentration is again somewhat 
on the high side (see Fig. [SI right panel). 

We have also checked that outside 100" radial 
anisotropy is disfavoured: even when forcing the T„ to lower 
values (we tried different values down to T» = 5), in or- 
der to allow for more radial anisotropy, the match to the 
outer data, especially the kurtosis, was poor (see dashed 
orange line). This result is somewhat surprising since pre- 
dictions from galaxy formation simulations generally show 
a significant degree of radial anisotropy (see e.g. ML05 and 
references therein), which has been confirmed by dynam- 
ical analysis in t he case of a few galaxies ( R+03; N+09; 
DL+09; but see iForestell fc Gebhardj lioTol ). Indeed, we 
have used directly the ML05 expression (see Eq. in mod- 
elling our data and found that the fit to both the VD and 
the kurtosis was possible only with a too small r a (~ 6"), 
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Figure 8. As Fig. [7] Confidence level of the p s — r s parameters marginalized with respect to Y, and /3o or r a (except for "NFW+ high 
/Jo" which is not a best-fit model). The "NFW+/3(r)" model is repeated as overlap with Fig. [7]. 



which is complet e ly inc onsistent with the values found by 
iMamon fc Lokasl |2005l . i.e. 1.4 R c ; see Fig. [8] gray dot- 
dashed line). Fixing r a to the expected value, the fit was 
possible only with a larger T» ~ 6.5. In either case, though, 
a much poorer significance of the fit than the one given by 
our preferred /3(r) profile (Eq. [5| was found. 

In summary, our exploration of the NFW models indi- 
cates that halo parameters corresponding to WMAP5 ex- 
pectations are compatible with the data. The agreement is 
better for a Salpeter IMF, with the concentration becoming 
somewhat high for a Kroupa IMF. The near-isotropic or- 
bital distribution that we infer is at odds with standard pre- 
dictions for radial orbits. However, as discussed in the Ap- 
pendix, there are some uncertainties in the classification of 
velocity outliers, such that we cannot yet claim the isotropy 
conclusion as robust. 

3.5.3 Effect of adiabatic contraction 

The baryonic collapse occurring during galaxy assembly is 
one of the physical process that can shape the central DM 
distribution in a way different from the predictions of the 
dark matter only N-body simulations. Given a dark matter 
halo distribution with the properties predicted by such simu- 
lations, the (collisional) collapsing gas can exert a dynamical 
drag on the DM particles and produce a more concentrated 
final DM density profile (see e.g. Blumenthal et al. 1986). 
The net effect is a larger central DM fraction and conse- 



quently a lower stellar mass contribution (i.e. a lower T„) 
to the total mass in the central regions (for fixed dynamical 
M/L and halo parameters). 

This process can be described analytically by an adi- 
abatic contraction (AC hereafter; iBlumenthal et all 1 19861 : 
iGnedin et al]|2004 G+04 hereafter) of the dark halo. Since 
there is not yet a final consensus on the effectiveness 
and accuracy of the descriptions on the market (see e.g . 
iPedrosa et "al]|2010l : [Duffy et alj|20ld : iTissera et aUl20ld ). 
we decided to use the recipe from G+04. The G+04 model 
produces a weaker effect on the final DM distribution than 
the original Blumenthal recipe, and appears closer to the re- 
sults obtained in the cosmological simulations including the 
baryon physics. 

A critical evaluation of the baryonic processes is beyond 
the purpose of this analysis, where we only intend to check 
whether including an analytical recipe for AC in our Jeans 
analysis would provide a viable explanation to reconcile the 
estimated T* derived from our analysis and the stellar pop- 
ulation models. Furthermore, to our knowledge, the use of 
the AC in detailed Jeans modelling of the velocity dispersion 
profile of an elliptical galaxy has not been attempted before, 
so we consider this an interesting exercise even though the 
AC recipe might be not optimal. 

For this purpose, in our equations the total mass gen- 
erating the potential $ = GM(r)/r is given by considering 
as an adiabatic invariant the quantity 
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Figure 9. As Figs. [jj and [8] with LOG models. The right-hand panel shows the corresponding 1 and 2 a confidence level of the vo — ro 
parameters marginalized with respect to T* and /3 parameters (when available). The curves correspond to models as in the panel legends. 
See text for details. 



M(f)r = const (18) 

where x = Ax™ and x — r/r v i r . By calibrating Eq. 1181 to 
collisional N-body simulations, G+04 have fixed A = 0.85 
and w = 0.8. The contracted DM mass distribution has been 
derived by solving the equation 

[M tot (f)]r = [M DM (f) + M t {r f )]r f (19) 

where M to t = Mdm + M*, and Mum and M* are the fi- 
nal dark and stellar mass respectively (initially assumed 
to have the same spatial distribution). The model results 
are shown in Fig. [5] and the model parameters in Table 
[2] ["NFW+AC+iso, +/3 , +0(r)"]. There are two main re- 
marks that we can derive from these results. 

First, since the effect of the AC is to drag more DM 
into the central regions, the Y„ turns out to be smaller 
than in the no-AC case. For the isotropic case we obtain 
T* = 5.7 Y@ t v (see Fig. [JJ tick purple dashed line), but if 
we again include /3(r) as in Eq. [6] the best fit is found for 
T* = 5.5 Tq^v and r a = 33". The goodness of these fits is 
slightly worse than, but similar to, the uncontracted NFW 
models (see Table [2]) with the model curves looking very 
similar to the eye (see Fig. [8] tick gray linej^l. 

Second, the (pre-contraction) dark halo parameters 

5 The model with constant anisotropy and AC yielded a rela- 
tively poor fit, and a very high halo concentration (Table [2{. 



turn out to be in very good agreement with ACDM. E.g., 
for the anisotropic model, the NFW dark halo turns out 
to have c v ir = 7.5 which matches the WMAP5 expectation 
(cwmaps ~ 7 for logAf v i r = 13.4). 

When forcing the fit to a lower T, = 5 Tq^v", the halo 
parameters change slightly: the best fit is c v i r ~ 9.1 and 
logAfvir = 13.5, which is higher than the typical prediction 
but still consistent with the scatter. 

This is one of the most notable results of this paper: 
for the first time using stellar kinematics extended out to 
~ 5 R e , it has been demonstrated that the dark matter 
content of a giant elliptical galaxy may be compatible with 
ACDM. 

3.6 Multi-component model results: LOG model 

We next carry out the model sequence for the LOG mass 
model (Section l3.2.2l) . with results shown in Fig. [9] and Ta- 
bled 

3.6.1 The isotropic model 

For the isotropic case, the LOG model can fit the data better 
than the NFW model in the central regions and equivalently 
well in the outer regions (see Fig.0 thin green dashed line). 
This is because the LOG potential has an internal core with 
little DM contribution in the central regions. In this case we 
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Figure 10. Radial mass distribution of NGC 4374. The left panel shows the cumulative mass, and the right panel shows th e circular 
velocity profile. Model curves from this work are as in the legend. We also show the v c profile from iKronawitter et al.l i2000h (shaded 
area includes the variance of their models) . 



find also a large stellar mass-to-light ratio, T* = 6.6 Tq,v, 
which is more compatible with a Salpeter IMF than Kroupa. 
A massive DM halo is required outside ~ 100" (mean vrj ~ 
450km s - ; see Table[2] "LOG iso"), consistently with the 
pseudo-mass inversion analysis and the NFW solution (see 

Fig. Eg. 



3.6.2 Models with orbital anisotropy 

Adopting a constant non-zero anisotropy (/3rj = 0.3) allows 
for a Kroupa-compatible T* = 5.5 T©,y (the same as found 
using the NFW+AC model). However, the fit is poorer (see 
Table[2] "LOG+/3o")i m particular at very small radii (even 
though these are penalized in our model) and owing to the 
higher estimates of the kurtosis at R > 100" (=2 dex), as 
shown in Fig. [9] (thin red dot-dashed line). 

We have checked if larger /3rj could be consistent with 
the data at large radii and found that once M/L* and r c are 
fixed, there is a degeneracy between the v c and the /3o values: 
a reasonable fit to the data is obtained for v c — 410 km s _1 
and /3rj = 0.1 and v c — 470 km s _1 and /3o = 0.5 with 
M/L* = 6 Tq^v and r c = 25 kpc. Once again, the kurto- 
sis helps to put constraints on the allowed /3o- the xVd.o.f. 
calculated over only the model versus observed kurtosis pro- 
files is much smaller for /3rj = 0.1 (~ 9/20) than for fio = 0.3 
(~ 12/20) and /3 = 0.5 (~ 22/20), which is a final demon- 
stration that strong anisotropy can be excluded at large 
radii. 

Finally, we have adopted the /3(r) as in Eq.[S] The best- 
fit model is not showed (but almost identical to the one with 
AC as in i]3.6.3[) . while parameters are reported in Table [2] 
["LOG+/3(r)"]. The anisotropy radius r a is slightly larger 
that the one estimated with the pseudo-inversion procedure 
and NFW (r n = 45"), although the /3(r) profile turns out 
to be almost unaltered. The Y* = 6Y0,y is closer to the 
isotropic case, since this is mainly constrained by the central 
regions which are almost isotropic according to the Eq. [6] 



3.6.3 Adiabatic Contraction 

For completeness, we have modelled the effects of a hypo- 
thetical AC on the LOG halo. Because of the non-cuspy 
nature of the initial halo, AC turns out to have only a weak 
affect, and does not change any of the above conclusions. 
Model curves are almost indistinguishable from the ones 
with no-AC as shown in Fig. [§] (green thick dashed line: 
isotropy; red thick dot-dashed line: constant anisotropy) as 
a consequence of best-fit parameters very close to the ones 
obtained for no-AC [Tabled "LOG+AC+iso, +/3 , +/3(r)", 
and confidence contours in Fig. [5] . 

Finally, the simultaneous use of the /3(r) anisotropy 
as in Eq. [6] and the AC allowed the best fit to the data 
(black thick line) as for the NFW case. For the LOG po- 
tential the stellar M/L turned out to be T, = 5.5T0,v and 
vo — 403 km s _1 (see Table [2j, and the anisotropy radius 
turned out to be very similar to the NFW models (r a = 35"). 
Once again the AC seemed to be a crucial ingredient to al- 
leviate the problem of the stellar M/L problem by naturally 
accommodating a Kroupa- like T„. 



3.7 Summarizing the best halo models: mass 
profiles and circular velocities 

Before we discuss the implications of the best-fit solutions 
from the previous sections, we summarize the models which 
we consider more physically meaningful. As seen in Table [2] 
and discussed earlier, most of the models presented are sta- 
tistically good fits (e.g. the reduced x 2 = xVd.o.f. is almost 
everywhere < 1), but some of the models were incompatible 
with related theoretical predictions. 

For instance, the no-AC models "NFW+/?o" and 
"NFW+j3(r)" have x 2 = 0.5, 0.4 respectively, but the im- 
plied halo concentrations are improbable given the ACDM 
expectations. Also "NFW+iso" has a rather small x 2 = 
0.6 and a fairly ACDM-like halo, but the large Salpeter- 
like Tq^ makes this solution unfavourable. On the other 
hand the model "NFW+AC+/3(r)" has a X 2 = 0.45 and 
is fully consistent with both ACDM concentrations and 
a Kroupa IMF, and so is considered our best reference 
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model. For similar reasons, the favoured LOG models are 
the "LOG+AC+/3(r)'\ "LOG+/3 " and "LOG+AC+/3 "— 
all having \ ~ 0.65 and a Tq^v compatible with a Kroupa 
IMF. 

Going to the comparison among the different potentials 
compatible with the stellar kinematics, in Fig. [TU] we plot 
the mass profiles of some of these model solutions in order 
to gain a general sense of the different halo solutions accom- 
modated by the data. 

Considering the mass profiles, M(r), for the different 
models discussed above, the DM halo models (NFW and 
LOG) are very different from the no-DM case, with the v c 
remaining much flatter with radius than the stellar model. 

The mass profile at 5 R c (~ 30 kpc) is remarkably 
similar for the NFW and LOG models, demonstrating that 
this quantity is well constrained by the data, independently 
of the details of the mass models. 

Despite the uncertainties, for the NFW case the mass 
profiles as well as the w c jr profiles differ in the very central 
regions when comparing the un-contracted solutions and the 
contracted haloes. The relative normalization between the 
stellar and halo masses changes due to the higher dark mass 
allowed by the AC for a given halo concentration before the 
contraction. For the LOG model, T* seems to be more de- 
generate with the /3 value in the central regions (in the sense 
that higher f} would allow smaller T*, see H3.2.2|) . Overall, 
the Vdic profiles (Tig. llO[) turn out to be fairly similar among 
the different models up to the last datapoint (~ 340"), and 
beyond, if the profiles are extrapolated more deeply into the 
halo regions. Furthermore, the mass profiles are remarkably 
similar to the results of the pseudo-inversion method (see 
Fig. E| . 

Finally, in Fig. [10] we compare our results with the v c irc 
profile from K+00, which is based on long-slit data extend- 
ing out to ~ 70". Focusing on our LOG +/3(r) solution which 
is the most equivalent to theirs, our results are identical in 
the very central regions, with a slight discrepancy at larger 
radii. Not e that the f c j r c from K+00 extrapolated to 300" 
(Fig.17 in iGerhard et all l200lh is significantly lower than 
our new profile based on more extended data and models. 

The asymptotic run of all the model curves in Fig. [TO] 
is remarkably tight which means that at intermediate scales 
(of the order of the r s scale of the NFW haloes) the overall 
galaxy mass is quite well constrained and the scatter in- 
troduced by the halo models and the allowed anisotropy is 
small. However, an important cross-check would be to ver- 
ify how these models might differ around the virial radius, 
where the NFW and LOG profiles are expected to differ sig- 
nificantly (although the extrapolated M v ir values in Table [2] 
do not differ much). 



4 DISCUSSION 

The dynamical solutions for the bright elliptical NGC 
4374 all clearly indicate that this galaxy is surrounded 
by a massive DM halo. DM haloes were also found in 
four ordinary ETGs studied using PNe (not all of these 
studies using PN con straints): NGC 3379 (R+03; DL+09; 
IWeiimans etafl l2009h: NGC 4697 ( DL+0 8); NGC 4494 
fN+09; iRodionov fc Athanassoulal 120101); NG C 821 
I Weiimans et al.ll2009l ; iForestell fc GebhardtJ |20ich . Apart 



from alternative gravity theories (e.g. iTiret et al.l 120071 ) . it 
seems clear that elliptical galaxies in general contain DM, 
and the question is how the DM profiles compare in detail 
to predictions. 

Some of the galaxies above were modelled with NFW 
haloes and some with LOG haloes, while NGC 4374 is the 
first of these cases where both were tried. Unfortunately, we 
were not able to discriminate between the two models, given 
the limitations of the simple Jeans approach which cannot 
fit the observations in great detail and requires somewhat ar- 
bitrary weighting of the data points. Interestingly, the two 
models do seem to prefer different T* values, corresponding 
to Kroupa and Salpeter IMFs for the NFW and LOG haloes, 
respectively. Adopting a prior on the IMF may then pro- 
vide more information about the DM profile, and vice- versa. 
More detailed modelling may also be able to discriminate be- 
tween these haloes on the basis of dynamics alone: e.g. with 
much less extensive data in a sample o f ETG s, but using 
Schwarzschild modelling, I Thomas et al l (|2007i ) found some 
suggestions that LOG haloes were preferred over NFW. 

Adopting the NFW halo model for the time being, it 
is important to test the inferred halo parameters (density 
and scale-length, or virial mass and concentration) against 
predictions from cosmological simulations. N+09 assembled 
the PN-based results as well as a heterogeneous sample of 
other mass results from the literature. We reproduce this 
mass-concentration plot in Fig. 1111 with the theoretical pre- 
diction updated for the WMAP5 cosmological parameters. 

Although the mass profile uncertainties for any single 
galaxy are too large to make definitive statements, when 
considering a handful of galaxies together, a remarkable pat- 
tern begins to emerge. The fast-rotator or ordinary ETGs 
(along with spiral galaxies) appear to have low concentration 
haloes, and the slow rotators to have high concentrations, 
with a possible zone of avoidance in between, corresponding 
to the theoretical predictions. With the shift to WMAP5 
predictions, the low concentrations become less of a prob- 
lem, and the high concentrations more. 

Adding NGC 4374 to the diagram confirms this picture 
with a PN-based slow rotator analysis for the first time. The 
NFW solution with a standard Kroupa IMF coincides with 
the high-concentration region previously found for slow ro- 
tators using somewhat similar analyses. However, the story 
changes with certain modifications to the models. If the IMF 
is forced to Salpeter, less central DM is permitted and the 
implied concentration decreases. Alternatively, the high cen- 
tral DM content could be due to AC, with the "original" 
concentration much lower, as illustrated by the modelling. 
In either of these cases, the halo concentration becomes con- 
sistent with ACDM predictions. 

Selecting a "heavy" IMF or including AC may thus 
generally solve the concentration crisis for slow rotators — 
but what about the fast rotators? Although we have not 
explicitly modelled these galaxies with AC, some general 
trends may be gleaned from the ACDM-bas ed toy models of 
iNapolitano. Romanowskv fc Tortoral j2010). Their Figs. 6 
and 11 illustrate that for ETGs of all masses, AC is ex- 
pected to dramatically increase the fraction of DM found in 
the central R c . This implies that if AC were included in the 
models of the fast rotators of Fig. 1111 the halo concentra- 
tions which are currently on the margin of consistency with 
theory would become problematically low. 
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Figure 11. Dark matter halo virial mass and concentration pa- 
rameters. Several reference solutions for NGC 4374 (large filled 
circles) are plotted as well as other data taken from N+09. The 
blue and gray curves with surrounding shaded regions are the 
WMAP5 and WMAP1 predictions, respectively. The green small 
dot with error bars is the "NFW iso" solution (see Table [2] the 
stellar M/L is consistent with a Salpeter IMF), the black small 
dot is "NFW+/3(r)" (corresponding to a Kroupa IMF), and the 
big black dot to our favoured model "NFW+AC+/3(r)" . From 
N+09: Triangles and boxes mark fast-rotator and slow-rotator 
ETGs, respectively. The small filled symbols mark detailed ETG 
dynamical results using PNe and GCs (including error bars, where 
available). The open symbols show the dynamics-based ETG re- 
sults from N+05, with error bars in the upper right corner showing 
the typical uncertainties. The dashed line shows the mean result 
for X-ray bright groups and clusters, the dot-dashed curve is an 
inference for late-type galaxies, and the dotted curve is the trend 
from weak lensing of all types of galaxies and groups (see N+09 
for details). 



An alternative scenario might be to adjust the IMFs of 
the fast rotators to be lighter than Kroupa (Salpeter is in- 
cidentally too heavy in general for this class of galaxies; cf. 
C+06). This would allow for more central DM and conceiv- 
ably increase the inferred concentrations. 

In order to move all the "observed" ETG halo concen- 
trations into reasonable agreement with the predictions, we 
arrive at the tentative solution that (1) the slow rotators 
have Salpeter IMFs or AC, and (2) the fast rotators have 
ultra-light IMFs or no AC. If (1) and (2) are fulfilled, then 
there may still be a systematic concentration offset between 
fast and slow rotators, but this would be small enough to be 
plausibly explained by differing collapse redshifts. 

This solution would present the very interesting possi- 
bility that the fast and slow rotators are dramatically differ- 
ent in either their IMFs or their halo contraction histories. 
Systematic transitions in these properties have been sug- 
gested for various reasons in the past, but they appear to go 
in the wrong direction. In the modern "downsizing" picture 




Figure 12. Cumulative dark matter fraction as a function of 
radius. Results for different observed and simulated galaxies are 
indicated with different colours and linestyles as in the legend, 
with the results for NGC 4374 in black. The errorbar marks the 
typical error for the dark matter fraction of NGC 4374 at 5 -Re. 



of galaxy formation (e.g. 


Nelan et al.l 20051; 


Thomas et al. 


2005;ICimatti et al.ll2006; 


Pannella et alj 2006 


; Graves et al. 



120071 : ICalura et al] 120081 ) , the more massive galaxies like 
NGC 4374 would have on average formed their stars earlier 
and more rapid ly than in the more ordinary elliptic a ls. As 
summarized by iNapolitano. Romanowskv fc Tortoral (|2010l . 
Sec. 4.4), the IMF in these conditions is thought to have 
been if anything lighter rather than heavier. 

Also, as summarized in N+09, it is thought that AC 
could be counteracted by rapid, clumpy and starbursty as- 
sembly histories, whil e AC classically implies s mooth, slow 
gaseous infall (see also lLackner fc Ostrikerl2010l ). These con- 
ditions would suggest that the spirals should have stronger 
AC, and galaxies li ke NGC 4374 should hav e weaker AC (a 
point also made bv lChen fc McGaughll2008T ). 

Returning to a less model-dependent view of the mass 
profiles, we plot the cumulative DM fraction versus radius 
in Fig. 1121 as also done in N+09. The model inferences 
for NGC 4374 as well as some of the ordinary ETGs are 
plotted, along with examples from galaxy formation simula- 
tions (bo th in a full cosmo l ogical context and using ad-hoc 
merg ers; iDekel et alj 120051 ; iNaab et al.ll2007l ; lOiiorbe et al.l 
120071 ). Drawing attention to the 5 Re reference value, we 
see that the DM fraction for NGC 4374 of ~ 0.7-0.8 is sig- 
nificantly larger than what was found so far for ordinary 
ellipticals (~ 0.4-0.5), and similar to what has been found 
for group- and cluster-central ellipti cals (~ 0.8-0.9 u sing X- 
ray rather than dynamical methods; iDas et al]|2010l ). These 
results bracket the simulations values of ~ 0.5-0.6. 

The DM fraction results within 1 R c in Fig.[l2]based on 
detailed dynamical modelling at first glance do not seem to 
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square with other recent results from the literature. Various 
combinations of dynamical, strong gravitational lensing, and 
stellar populations analyses have found typical DM fractions 
within 1 R e of ~ 0.4 for fainter ellipticals and ~ 0.6 for 



20031; iTortora et al.ll 2009l. size and surface brightness (e.g 



Capaccioli et al. 19921). luminosity and velocity dispersion 
19761) a nd the colour/population proper- 



I Fabcr & Jackson 



brighter ones ([Napolitano. Romanowskv fc Tortora 20ld : 



ISchulz et aLll2010l ; lTortora et al.ll2010bl ; lAuger et alj|2010h . 
versus ~ 0.05 and ~ 0.3 here. 

However, in the case of NGC 4374, the ambiguity in 
the R e comes into play. In NRT10, galaxies of the same 
stellar mass have R c ~ 12 kpc on average, or ~ 145" at 
the distance of NGC 4374. Using this R c scale, we would 
have a DM fraction of ~ 0.5, consistent with the literature. 
As for the lower-luminosity ellipticals, NRT10 did find a 
fraction of galaxies (particularly ones with older stars) to 
have DM fractions lower than ~ 0.1, so the critical goal is to 
assemble a large sample of galaxies with detailed dynamical 
m odels to establish the trend s with good statistics. Fig. 1 
of lTruiillo -Gomez et al. (2010) does suggest that these three 
galaxies may happen to represent one extreme from a broad 
distribution of DM properties at intermediate luminosities. 
If this situation is true, the arguments above about halo 
concentration offsets would no longer apply. 



5 CONCLUSIONS 

We have presented a full Jeans analysis of the bright, slow- 
rotator elliptical NGC 4374 based on the observations of 
~ 450 PNe with the Planetary Nebula Spectrograph. The 
PN line-of-sight velocities extend out to ~ 5R e - We have 
constructed spherical Jeans dynamical models of the sys- 
tem: a "pseudo-inversion" model and multi-component mass 
models with and fourth-order moments constraints on the 
orbital anisotropy. 

The two approaches return similar values of M/L and 
anisotropy (see Fig. [5] and Table [5J) and both imply that 
NGC 4374 is a very dark matter dominated system with 
a near-isotropic orbital distribution in its halo. Dynamical 
analyses of more ordinary ETGs have suggested radially- 
biased anisotropy in their haloes as predicted by simula- 
tions (see Section [3X2J- The NGC 4374 result on the other 
hand would build on previous suggestions that slow rota- 
tors have surprisingly isotropic haloes, which would sug- 
gest a new scenario for building of the ext ended stellar 
envelopes of these galaxies ma y be required (jHwang et al.l 
|200S| ; iRomanowskv et al.l 1200911 . However, in this case the 
anisotropy result is sensitive to the assumptions about out- 
lier velocities, and further investigation is required. 

The mass profile results are on the other hand fairly 
insensitive to the outliers. The high DM fraction inferred 
within ~ 5 R c confirms the apparent dichotomy in DM 
content between slow and f ast rota tors proposed by N+09 
(see also iBertin et al.lll994l ; C+06; iNapolitano et alj|200j; 
C+09), and yields two important implications: (1) the DM 
dichotomy is not a result of systematic differences in the 
mass tracers used; (2) it is not a simple difference of group- 
central versus satellite galaxies since NGC 4374 does not 
appear to be at a group center (while the low-DM system 
NGC 3379 is). 

This apparent DM bimodality may mirror other transi- 
tions in ETG properties at similar luminosity scales, such 
as the relations between size and mass (e.g. IShen et al.l 



ties ( Tortora et al 



2010a) 



Given the limitations of the Jeans models and the stel- 
lar/dark mass degeneracy, we are not able to distinguish 
between different DM radial profiles, including LOG, NFW 
and NFW+AC haloes. The LOG models prefer high stellar 
masses consistent with a Salpeter IMF, NFW works with 
either Salpeter or Kroupa, and NFW+AC requires Kroupa. 
The nominal NFW+Kroupa model implies a halo with a 
concentration that is somewhat high, given WMAP5-based 
predictions. Adopting either Salpeter IMF or AC brings the 
inferred concentration down to the expected value. Thus, 
considering that AC has commonly been considered the de- 
fault expectation in galaxy formation, we have finally found 
an ETG analyzed using PNe that is naturally consistent 
with theoretical expectations for the DM halo. 

Comparing the NFW halo parameters obtained for 
NGC 4374 as well as for an assortment of other galaxies 
in the literature, we find evidence for the slow rotators to 
have much higher halo concentrations on average than the 
fast rotators. We discuss some possible variations in IMF 
and AC which could explain this difference, but there are 
also suggestions that the sample of fast rotator galaxies is a 
statistical fluke. 

Two primary avenues are needed to make further head- 
way in pinning down the properties of DM haloes in ETGs. 
One is to carry out more detailed dynamical and stellar pop- 
ulations modelling in an attempt to discern the DM profiles 
in detail. The other is to expand the sample of galaxies stud- 
ied, particularly at intermediate luminosities (Mb ~ —20). 
Work on both fronts is underway as part of the PN.S Ellip- 
tical Galaxy Survey. 
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APPENDIX A: ALTERNATIVE OUTLIER 
SELECTION AND DYNAMICAL 
IMPLICATIONS 

As discussed at the beginning of Section [2] a handful of 
"outlier" PNe were rejected from the overall sample using a 
3 <r "friendless" analysis. Although relatively few in number, 
the inclusion or exclusion of these objects in our dynamical 
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analyses could have a large impact on the conclusions, which 
we consider here in more detail. 

Fig. Q] showed the six outliers identified through this 
process. Two of them are extreme outliers and can be se- 
curely rejected, but the other four are only barely excluded 
at 3 a. This is a concern since in a data set of 450 ob- 
jects with a Gaussian velocity distribution, there should on 
average be one random object found past 3 a, and if a non- 
Gaussian distribution is allowed, then many more would be 
possiblQ 

We are not at a complete impasse since we notice that 
all the outliers have negative velocities relative to NGC 4374, 
which is not likely to be just a chance occurrenc 

We look at the situation in two-dimensions in Fig. IA1I 
focusing on the most outlying velocities. It turns out that 
the four most extreme velocities lie on an axis to the East of 
the galaxy's center, which is also the direction of the nearby 
giant elliptical NGC 4406 (M86) found ~ 1000" away. A 
similar pattern has been found in the globular cluster system 
of NGC 4374 (B. Kumar et al, in prep). 

The —1300 km s _1 relative systemic velocity of 
NGC 4406 (ti re i) provides a handy explanation for the 
low-velocity outliers — whether these objects simply belong 
to the NGC 4406 PN system seen in projection, or are 
part of an interactio n region between the two galaxies 
l|Arnaboldi et alj|l996h . 

Although a full analysis for such an interaction scenario 
is outside the scope of this paper, we can quantify the effect 
of a fly-by encounter between the two galaxies using the im- 
pulse approximation to estimate the energy injection int o the 
outer galaxy envelope (see e.g. iNapolitano et ai1l2002t ). We 
calculate an upper limit to this energy by assuming a tan- 
gential encounter with an impact parameter of b = 1000" : 



AE = 



4GMi_Mf- 



(Al) 



where G is the gravitational constant, M\ ~ 6 x 10 1 Mq 
(e.g., from the "NFW+iso2" model) is the mass of the per- 
turbed system (NGC 4374) which has been calculated within 
the impact parameter b, M2 = 0.5 x Mi at the same radius, 
and the mean square radius of NGC 4374, r 2 is taken as 
equivalent to the square of characteristic scale of the dark 
matter halo (~ 6.4 x 10 3 kpc 2 ). 

The resulting energy change is AE = 9.2 x 10 16 M 
km 2 s~ 2 which provides a heating contribution to the dis- 
persion of Cheat = 2A_B/(3A/ S h c ii), where M s h c ii is the mass 
of the galaxy shell which has experienced the energy trans- 
fer. Taking this shell in the radial range of 200"-1000" (i.e. 
> 3 R e ), we find (Xheat ~ 100 km s _1 . 

This extra heating term could handily explain the 
higher dispersion on the low velocity side as implied by the 



16 We have checked how the outlier velocities compare to 
the local escape velocities in our best-fit NFW halo (e.g. 
the NFW+AC+/3(r) in Table E}, which turns out to be 
~ 1250 km s — 1 relative to the systemic velocity. The two most 
extreme outliers would in this case not be bound to NGC 4374, 
but the other four could be. 

17 This asymmetry does not appear to be caused by an error in 
the adopted value of v ByB , as the peak of the LOSVD coincides 
with our self-consistent v ByB , which is in turn very close to the 
NED value. 
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Figure Al. Diagram of potential velocity outlier PNe. The 2D 
positions on the sky are shown relative to the center of NGC 4374, 
with the stellar isophotes at R m = 197" and 257" show as dashed 
ellipses. Squares represent approaching velocities and crosses are 
receding velocities, with symbol size proportional to relative ve- 
locity amplitude. The ~ 30 most extreme velocity PNe are shown 
along with the candidate outliers, to illustrate the normal velocity 
field of NGC 4374. The remaining PNe are shown as small grey 
points. The North and East directions are shown in the top-left 
corner. 



"outliers" in Fig. \T\ In this scenario, the close passage be- 
tween the galaxies would have heated the Eastern side of 
NGC 4374, with this event happening less than one-crossing 
time ago so that the asymmetry is preserved. Removing the 
four "outliers" would then restore the observed kinematics 
of the system to the approximate pre-interaction state, suit- 
able for equilibrium dynamical analyses. 

The interaction calculation above has been done un- 
der the assumption of the closest encounter (and the high- 
est energetic) allowed by the observed geometry. Any other 
less favourable configuration would produce a smaller en- 
ergy transfer and a more local effect of the encounter. In 
this case, the four Eastern low-velocity objects are likely to 
be true outliers, and the remaining two outliers to the North 
are less certain, and could be part of the normal velocity dis- 
tribution of NGC 4374. 

If those two objects are kept in the final sample then 
the velocity dispersion and kurtosis profiles are somewhat 
changed in the outer regions, as shown in Fig. IA2I The dis- 
persion profile becomes slightly flatter (slope —0.03 ± 0.07 
instead of the —0.07 found in £|2. 2|) . The kurtosis profile rises 
at large radii, where if we were to again use the equation B10 
approximation from N+09, we would infer a higher radial 
anisotropy (/3 ~ +0.4 instead of ~ —0.1). 

Carrying out some dynamical models as in the main 
Sections, we show in Fig. lA2l the results for the isotropic and 
j3{r) NFW mass models. We find best-fitting halo parame- 
ters of p s = 0.0019, 0.0030 M @ pc" 3 and r s = 110, 87 kpc re- 
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spectively, corresponding to c v i r = 8jl|, 101® an d l°g M v i T ~ 
13.6M0. These parameters are very similar to those found 
using our default outlier selection (see Figs. [7|and[8|, al- 
though the \ 2 fits are poorer. 

We also try out a more strongly varying /3(r) function 
motivated by the higher kurtosis, with high radial anisotropy 
at larger radii as illustrated by the right-hand panel of 
Fig. IA2I and named /J3 and /?4 some fixed profiles which 
bracket the tentative anisotropy value in the latest radial 
bin estimated as above (~ +0.4). In this case the fit is per- 
formed on the dispersion curve only. 

The quality of the corresponding dynamical model fit 
(top-left panel) is similar to the previous case, but the best- 
fit dark matter halo turns out to be almost identical for the 
two anisotropy profiles and have a higher halo concentra- 
tion and small virial mass, slightly off the WMAP5-ACDM 
predictions: p s = 0.006 Mq pc -3 and r s — 46 kpc, cor- 
responding to c v i r = 13lg and logM v i r ~ 13.1 ± O.IM©, 
whith errors including the variance of the assumed /3 pro- 
files. In both these strongly radial models, the velocity dis- 
persion bends quite significantly outside the last dispersion 
bin, which is a prediction that should be tested with more 
extended data. 

We thus find that the impact of the outlier ambiguity 
is confined to the anisotropy conclusions, with highly radial 
halo orbits suggested by the kurtosis but hardly matched by 
the dispersion profile which is flatter when the two uncertain 
outliers are included. The mass profile inferences are pre- 
sumably unaffected because of the pinch-point phenomenon, 
whereby the projected dispersion is only weakly dependent 
on anisotropy in certain regions of the gal axy. Further ob- 
serva tions of PNe at larger radii (see e.g. lArnaboldi et al.l 
particularly on the West side of the galaxy, could 
clarify the situation by more strongly constraining the dis- 
persion and kurtosis profiles past the pinch point. 
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Figure A2. Effect of the outliers selection. Top: the velocity dispersion profile of the outliers selection with the friendless algorithm, 
adopted in the paper (black) is compared with the dispersion obtained by including the two uncertain outliers along the 3cr borderline 
(red). Bottom: the same for the kurtosis profile. Overplotted there are the isotropic model obtained for both profiles, and a best fit model 
with a steeply increasing (3(r) profile as suggested by the kurtosis including the two uncertain outliers, as in the top-right panel. Models 
are as in the legend. See text for details. 
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